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ON THE CHARACTERISTIC EXPONENTS OF FLOQUET 
SOLUTIONS TO THE MATHIEU EQUATION 

JAN ERIC STRANG 

Abstract. We study the Floquet solutions of the Mathieu equation. In order 
to find an explicit relation between the characteristic exponents and their 
corresponding eigenvalues of the Mathieu operator, we consider the Whittaker- 
Hill formula. This gives an explicit relation between the eigenvalue and its 
characteristic exponent. The equation is explicit up to a determinant of an 
infinite dimensional matrix. We find a third-order linear recursion for which 
this determinant is exactly the limit. An explicit solution for third-order linear 
recursions is obtained which enables us to write the determinant explicitly. 



1. Introduction 
The Mathieu equation is a special case of Hill's equation 

^ + (a-2#(z))^ = 

where 4> {z) is a 7r-periodic differentiable function of the real variable z of maximum 
absolute value 1 and a, q are real constant^. Thus, the equation can be written in 
terms of the Fourier transform of 4> (z) as 

(1.1) ^0 o + 2^02rCOs(2r.z)^ = O 

if the function <j> (z) is even. The Mathieu equation is a restriction to <j> ( z ) — cos (2z) 
of the Hill equation, i.e. 

(1.2) ^j + (a-2 9 cos(2z))V = 0. 

The equation was originally studied in the context of vibrational modes of an 
elliptic membrane by E. Mathieu pQ. Considering time periodic solutions, the 
initial wave equation can be separated and written in elliptic coordinates [3] , 
consequently the wave equation takes the form of eq. (|1.2p . In a similar fashion, this 
procedure can be applied to heat equations on elliptic domains [2]. Furthermore, 
the equation appears in solid mechanics with a periodic perturbation, e.g. for a rod 
fixed on one end and being under periodic tension at the other end [2]. 
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throughout this paper we will make an effort to use existing conventions regarding the nota- 
tions of the Hill and Mathieu equations (see for example [2]) 



The Hill equation also appears in quantum mechanics as an eigenvalue problem 
for a particle in a periodic potential [3], If we restrict the potential to 2Qcos(2z) 
in one dimension, the Schrodinger equation (h = 1, 2m = 1) 

dib d 2 ib n ^ . . , 
*W = "9# + 2QC ° S(2Z) ^ 
will lead to an eigenvalue problem for the operator H = —d 2 /dz 2 + 2Q cos (2z), 

£ty> = + 2Qcos(2z)i/', 

which is easily written in the form of eq. I|1.2p 

(1.3) ^| + (£-2Qcos(2z))</> = 0, 

identifying a = E and q = Q. This particular problem is of special importance in 
solid state physics where the crystal lattice atoms constitute a periodic potential 
for loosely bound electrons. Eq. l|1.3p is a crude one dimensional simplification of 
such a potential. 

After a short introduction to the Mathieu equation, we will describe (section 
12. ip the Floquet theorem and its application to the Mathieu equation. We will 
then state some general properties of the Floquet solutions (section I2.2|) and more 
specifically the Whittaker-Hill formula (section I2.3p which has given rise to this 
work. We will then focus on the calculation of an infinite dimensional determinant 
(section [3]) needed to solve the equation given by the Whittaker-Hill formula. In 
section [3721 we describe a method to obtain an explicit result for a third order linear 
recursion. In section l3~3l an explicit form of this determinant is given. 

2. The Mathieu equation 

There exist a multitude of solutions to the Mathieu equation (|1.2p . both periodic 
and non-periodic. What is distinctive is that they cannot generally be written 
in terms of elementary functions, i.e., they are transcendental functions. For a 
summary see, e.g., 0E] and for a more complete account with examples see [3]. 
There also exists a number of approximation schemes for the Floquet solutions jH 
O O [3]. To our knowledge, the most recent approximations are due to Frenkel and 
Portugal [7|. We will restrict our study to periodic solutions and more specifically 
to the Floquet solutions of eq. (|1.2p . 

2.1. Floquet Solutions. As the coefficients of the Mathieu eq. (|1.2p are periodic, 
one can apply the Floquet theorem. According to Floquet [8], a differential equation 
with periodic coefficients admits at least one periodic solution of the second kind, 
i.e., if tp (z) is a uniformly continuous solution, then 

(2.1) if) (z + r; a, q) = eip (z; a, q) , Vz G K, 

where r is the period of the coefficients and the multiplier e is a constant. If e = 1 
the solution is of the same period as the coefficients and is called periodic of the 
first kind. Further identification of the multiplier leads to the solutions of the form 

(2.2) V fa a, q) = e tl ' z u lJ (z; a, q) , v G E, 

where u v (z) is a periodic and uniformly continuous function of the first kind with 
the same period as the coefficient. We will call such solutions Floquet solutions to 
the Mathieu equation, and call v its characteristic exponent. 

It should be pointed out that though the sum of Floquet solutions is a solution to 
eq. I|1.2p . it is no longer a Floquet solution. Furthermore, the solutions according to 
eq. i|2.2p with real v constitute so-called stable solutions. There exists values (a, q) 
for which such solutions do not exist, in which case v becomes complex. These 
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solutions are no longer bounded and hence are called unstable. The eigenvalues a 
for which v is real form bands [9] , with gaps corresponding to complex characteristic 
exponents. For non integer characteristic exponent, there exists a space of solutions 
spanned by two independent solutions [H El GEE]- F° r integer values of v there exists 
only one solution which is usually written in terms of the Mathieu function me„ 
OElll, 2j. Integer v are border solutions to the eigenvalue bands. 

We will state an important connection between the eigenvalue and the charac- 
teristic exponent as given by Kohn [I]. The solutions ipi, ^2 to a given pair (a, q) 
with initial conditions 

i> x (0; a, q) = 1, $ (0; a, q) = 

^2 (0; a, q) = 0, ip' 2 (0; a,q) = l 

form a fundamental system, i.e., they are linearly independent [H|, Any solution to 
a pair (a, q) can be written with those solutions. A Floquet solution satisfies 

ip (z + t; a, q) = eijj (z; a, q) , 

hence 

ip(r;a,q) = eip(0;a,q), 
ip' (r; a, q) = Eip' (0; a, q) , 
where e = e WT . Rewriting ip (z; a, q) as a linear combination of ipi ( z ; a i q) an d 

-02 0; «, q) 

1p=ATp l + Bip 2 

one obtains the equation 

(2.3) e 2 - 2//(a)e + 1 = 

upon elimination of A and B. Here fi (a) is an entire function of a that can be 
expressed as 

m ( a ) = 2 ( r ' a ' + ^2 ( t ; a > ?)) • 

Insertion of e — e %vr according to the Floquet theorem in eq. (|2.3p and subsequent 
multiplication by e~ lUT gives 

(2.4) cos(j/t) = fi(a) . 

This equation can be used to obtain a relation between the characteristic exponents 
and the eigenvalues of the operator H. Unfortunately the function fi (a) cannot be 
written in such a way as to solve eq. (|2.4p . 

2.2. Determinantal equation. Since u u is periodic, we will seek solutions u v £ 
L 2 ([0, 7r]) and extend these solution to z S R. We may expand the entire solution 
as the Fourier series 

k£Z 

which, inserted into eq. (|1.2p . gives the recurrence relation 

(2.5) {{2k - vf - aj c 2k + q (c 2 («+i) + c 2 ( K -i)) = 0, Vk G Z 

for the Fourier coefficients C2 K (a, ^)- Since u v 6 L 2 ([0, 7r]), we see that {c2n} Ke z € 
^2 (^) whereupon 

yj |c2« (^; a, <?)| 2 < 00. 

For any finite (truncated) upper limit neZ, eq. I|2.5p can be written as a matrix 
equation, 

(2.6) A n (p; a, q) c„ = 0. 
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where c„ = (c_ 2 „, ■ ■ ■ , c_ 2 , c , c 2 , c 4 , . . . , c 2 „) and A„ (i/; a, q) is a 2n + l x 2n + l 
matrix. For the sake of convergence (see below) we will define the equivalent system 
of equations 



(2.7) 



C2k 



1 (c 2 ( K +l) + C 2 ( re _l)) 



(2k — v) — a 
whereupon eq. I|2.5p takes the form 



o, Vk e z, 



(2.8) 

with 
(2.9) 



An (u; a, q) 



A n (u; a, q) c„ = 



/ 1 

£,2n-2 



£4 


1 


£4 

















£2 


1 



















£o 


1 


£o 

















£-2 


1 


i-2 

















£-4 


1 


£-4 



-2/1, 



£-2ri+2 
1 



and 



^2k 



(2k- i/) 

For the finite dimensional case, finding the non-trivial solutions of eq. (|2.8p is 
equivalent to demanding 

det (An) = 0. 

The determinant can, however, not always be defined for the infinite dimensional 
case. 

The determinant det (C + I) of an operator C + 1 defined on an infinite dimen- 
sional Hilbert space, can in certain cases be defined by the product of (1 + Aj) 
where \ are the eigenvalues of C. For trace class operators C, this product is well 
defined, independent of the choice of basis and converges absolutely [10]. A trace 
class operator, suitably defined on a separable Hilbert space, has the property 



Tr(|C|)=^/Xi<oo ) 



where /Ltj = |Aj| are the singular values of C |10[ fTT] . For such operators one can 
define the trace norm 

oo 

ii c Htt =X^- 

We can consider the operator B n : C 2n+1 -> C 2n+1 with B n := A n - I . By 
letting n — ► oo , the limiting matrices .4, B will constitute operators on the infinite 
dimensional Hilbert space £2 (Z). Furthermore, *6 can be written as 

B = B +1 

where B+i and are single sub- and super-diagonal, respectively, with elements 
(• ■ ■ , £-2n, £-2n+2, ■ ■ ■ , £-2, £o, £2, ■ ■ ■)• The singular values of these operators will 
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hence be \%2k\ , k € Z, The sum of the singular values of these operators will be 



Tr(|B ±1 |) = £|6„| 

q 

(2k — v) 2 — a 



which can be rewritten 

(2-io) Ei^i = ~£ 

This sum will converge since Y1T «*» X > lj converges and the terms under the 
nominator in eq. (|2.10p are inferior to \r\~ p , 1 < p < 2 for large enough k. This is 
valid for all finite values of a and v for which no denominator of vanishes, i.e., 

(2.11) 2k-u^±VE. 

The operators B+\ and B-\ defined on £ 2 (Z) are therefore trace class operators 
with trace norms 

||B ±1 || lv = Tr(|B ±1 |). 
Furthermore, by the triangle inequality of the trace norm [TT] we can write 

<||B + i||av+l|B-illiv, 

thus 

\\B\\ Tr <2-\\B ±1 \\ Tr <cx>. 

The operator B is hence also of trace class. By the properties of trace class operators 
[T2l fTOl fTT] one can define the convergent determinant det (A) = Yl K ez A K (A) of 
the operator I + B = A where A K (A) are the eigenvalues of A. The operator, by 
the Fredholm alternatives, is then boundedly invertible if and only if det (A) i= 

E2|. 

The recursion system of eq. i|2.5p will render ip zero if two consecutive coefficients 
vanish. Seeking the non-trivial solutions of eq. (|2.6p is then equivalent to demanding 
that 

det (A(v; a, qj) = 0. 

2.3. Whittaker-Hill formula. Whittaker was able to produce an astounding for- 
mula for the determinant corresponding to the full Hill equation [14]. We will 
describe this result for the determinant A (a, v) = det (A [y; a, q)) for the Mathieu 
problem following the discourse of McLachlan [2] . 

Since k takes all values from —00 to +00, the determinant A (a, v) will be left 
invariant under the change v — > — v, hence A (a, v) is even in u, i.e., A (a, —v) = 
A (a, v). Furthermore, (2k-(i/ + 2)) = (2(k+1) - v), giving A (a, v + 1) 
A : r ; . v) again since k takes all values from —00 to +00, i.e. the determinant is 
a periodic function in v with period 1. We can hence restrict the study of the 
analytic behavior of A (a, —v) to the strip < v < 1. A is a determinant and thus 
is a sum of products of & (a, z/). Disregarding their poles, the functions & (a, z/) 
are analytic functions hence so is A except at the poles of the functions & (a, 1/). 
We see by the structure of *4 that no term in £2^ will appear in higher powers than 
1, i.e., all poles for the values of a and v satisfying eq. (|2.1ip will be simple. 

The function V (a, is), 



cos (ttv) — cos (ny/a) 
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also has simple poles at those values of a and v. If the function C (u) is appropriately 
chosen, the function 

9 (a, v) = A (a, v) - C (u) ■ V (a, v) 

has no singularities. Note that A (a, v) has only one pole on the strip < v < 1. By 
the periodicity (in v) of A (a, v), finding the function C (v) is reduced to calculating 
the constant C = C (v) corresponding to the quotient between the residuals of the 
functions A (a, v) and V (a, v) at that pole. With the proper value of C, 6 (a, v) is 
analytic on the whole complex plane and has no poles. It must hence, by Liouville's 
theorem, be a constant. 

In the limiting case v — > +ioo, the matrix A (v] a, q) becomes diagonal dominant, 
i.e., only the diagonal elements (equal to 1) remain in this limit since lim^^^ £ 2 k — 
0, Vk € Z. Thus 

lim A (a, v) = 1. 

We see further that Ym\ v ^, loa V (a, v) — 0, hence 6 (a, v) = \. It follows that 

A (a, v)-\ 

2-12 C= I ' \ ■ 

For the case z/ = we have I? (a, 0) = 1/ (1 — cos (Tr^/a)). Excluding the values of 
a = 4k 2 for which A (a, v = 0) cannot be defined, we can write 

c ^ A(q,0)-1 

1 — cos (ny/a) 

By substitution into eq. (|2.12p one finds that 

2 (™\ = A(n nUm 2 



(2.13) sin^ ( T j = A (a, 0) sin^ , a ^ 4 K 2 . 

The characteristic exponents resulting from the eigenvalue a are solutions of this 
equation. This formula was shown by Whittaker [2] for the full Hill problem. The 
full Hill problem of eq. Ijl.ip has an infinite number of Fourier coefficients. The 
matrix A will in this case be a full matrix, i.e., not merely a tridiagonal matrix as 
in the Mathieu problem. 

3. Calculating A (0) 

3.1. A (0) as the limiting case of a recursion. By the Whittaker-Hill formula, 
the calculation of the characteristic exponents of the eigenvalue a is reduced to 
the calculation of A (0) = A (a, v = 0). The determinant of an infinite dimensional 
square matrix can be defined by recursion. Defining the square tridiagonal matrices 

/ 1 & \ 



J2 = 



t0 i tO 

S2(i-1) 1 S2(<-1) 
CO 

»2(i-2) 



e5r; ■» 1 



1 $(i-2) 

tO i cO 

V o e 2l i / 

with diagonal elements 1. The off diagonal elements are collected in the vectors 

(gU, C°- 2 , ••••^ and (4°, ....CI £ (i _ 2) , e (< _i)). Here de- 

notes ^ (z/ = 0) as previously defined. The determinant A (0) is hence defined as 



the limiting case of A? = det L/4°) , i.e., 

A(0) = lim A°. 

i — >oo 

Since for j/ = 0, = we see that the matrices A% are symmetric with respect 
to the second diagonal. 

We see that any matrix A®, i > 2, can be written in terms of A®_ 1 through 

/ 1 g \ 



^4° 



V 



i > 2, 



$ i y 

where the dotted square represent the outer elements of A®_ 1 . Defining the trun- 
cation operators of a matrix M; r(M), l(M), u(M), d(M), i.e., the matrices M 
without left (r(M)), respectively right (l(M)), columns and the upper (d(M)), 
respectively lower (u (M)), rows, we can write the Laplace decomposition of the 
determinant of the matrix A®. Developing along the first row we get 



det (.4?) 



A , 



-c 



and then developing the first term on the last row and the second term by the first 
column, 



det (AT, 



det (AU) - e? 



Si Si-1 



dr(^_ x ) 



to 

Si-1 



Sz 



Si-1 
1 



Proceeding in a similar fashion one can write 

det(A°) = AO_ 1 -e i o e i det( u /(^_ 1 ))-^ 1 det(rd(^_ 1 ))- 
+ mU) 2 det (uldriAU)). 

One sees directly that uldr (A®_i) = -4°_ 2 . Furthermore, by the symmetry prop- 
erties around the second diagonal, det (rd (A®^)) = det (yd (A®_i)) . Thus we 
conclude that 

(3.1) A° = A?_! - 2 • det (rd (4°^)) + (^if A°_2- 
By simply rewriting eq. (|3.ip 

A° - A°_ x = (tf&J 2 A?_ 2 - 2 • det (rd ■ 

we see that the series converges since Hindoo £i = 0. 

By further defining £1; = det (yd (A®)) — det (rd (A®)) and decomposing by 
minors, one can write 

(3.2) n i = det(Al 1 )-^_ 1 n i - 1 . 



Substituting det (rd in eq. (|3.1|l one finds the recursion in A», 

(3.3) A 4 = /3iAi-i - ai/3jAj_ 2 + A<_3, 

where we have substituted = and 1 — a, = /3». 

3.2. A solution to linear third-order recursions without constant terms. 

We consider the recursion 



(3.4) 



By simply rewriting it for i — 1 and multiplicating it with —a/bi-i, assuming 



(3.5) 



Summing eqs. (|3.4|l and (|3.5p one finds 

c 



Cit^ X ^% X 

Ui-i = : Ui-2 - CiUi- 3 Ui-4. 



(3.6) 



ut-i + [bi - 



h-i 



(■iCi—l 

Ui-2 ; Ui-i. 



One can hence, by this iteration procedure, construct a sum equal to some Ui over 
n elements u&, where the coefficient for Ui_„ + i vanishes, 



(3.7) 



Uk-i 



{Qn}u k - 



Here {g„} k _- are the coefficients corresponding to the sum equal to u& with last 
term in the sum Uk—n- Using this ansatz and eliminating {G%} Uk—n as previously 



one finds that 



{9n+l} 



k-j 



{9n}k-j 

\ y n> k-(n-2) ~ 6*_(„_ 3 ) 



1 < i < n — 3, 



j = n-2 



{G k n+1 } = -K} 



&fe-(n-2) 



Since these recursions are of first order, one can deduce 
f 0*}, 



{ffn+l}fc-j - 



1 fc-j 

(-ir 5 -{^} fe _ ( „_ 2) 



X 



n— 6 "fc-( ra -2) ,c fc ,c fc-i 

&fc-(n-2)6fc-l 

n— 3 c fc _j 



(-1) 



, 1 < i < n - 3 
,j = n-2 

, j = n - 1 



x n j= 3 6^7 



{^n+l} 



n— 3 

r jyi— 5 c fc Cfc— i tt c k-j 
bk-i +\ b k -j 

3=3 J 



from the initial values of the recursion given by eq. (|3.6|) . We can hence write all 
{9n} k -j and {Sn+i} explicitly. 

Being able to write such sums for any relevant Uk and n, one can eliminate the 
term of highest index on the r.h.s. of eq. (|3.7|) . If we consider 



Uk = 



n-2 



+ {Gn}u k - 



3=2 



{Q k n}u k - 



and 



Uk-l 



E/ {fim }fe_l- 



U k - 



3 = 1 



summing the two we find 



Uk = 



n-2 

E {a^k-jUk-j 

3=2 



E {9m l } k -j U k-i 



3=2 



+ {s£} fc _ 1 {<&- 1 }«*-i- 

If we choose m such that 



find 



m + 1 = n, 

U k = ({<tf} fc _ 2 + {ff£} fc _ 1 {<fc 1 l} fc _ 2 )«fc-2 + 

n-2 

E {{9n}k-j + {9n)k-l iatDk-j) U ^-3 



3=3 



+ 



(3.8) 



+ {m + {9t} k - 1 K-\})uk-r, 



Upon t such eliminations, one can hence write Uk as 

n-2 

(3-9) E HYk-^k-i + iTlYuk-n 

j=e+i 

where {lnY k -j an< ^ {^«}^ are the coefficients corresponding to Uk-j and u k - n re- 
spectively, after i eliminations. If we proceed as previously and eliminate u k _^ +1 ^, 
we get 

n-2 

«* = E ({7"}L J + K}L-i{^} fe _,)^-i + 



Ufc- 



+ (W}L-i{«} + {r*} < 

Thus we can define recursions for {7*}^ . and {rjj} £ through 

(3-io) HY k + _) - mU + mlt-i > 



which are second, respectively first, order recurrences. 

By calculating the coefficient of u k -j of highest index, i.e., {in} k-e-i' one can 
recognize a structure for these coefficients. {lnYk-e-i can ^ e wr itten as a sum of 
products of < g k n ~Jt \ where the number of terms in the products range over all 

I J k — v 



values from 1 to I + 1. As an example, we consider such a sum for {7^}, 7 for 
terms with four <^ g^Z!t \ factors. These terms are written in the left column and 

I J k — v 

the indices (fi, v) corresponding to the terms I \ are shown in the right 
column, 





-1 \y-a- 


-n 
-i/fc- 


-2 Ifn- 


1) 

-2/ fe- 


{q k ~ 

-3 


11 + 

-3/ fc-6 T 






+ Ui k ) 


-1 \iin- 


-n 
-1/ 


-2 l^n- 


ll 
-2/fc- 


-4 \Hn- 


-4/fc-6 ^ 


: (0 1) 


fl 2) (2 4) (4 6) 


+ U k \ 


\q k ' 

-1 


-n 
-i/fe- 


-2 \ y n- 


-2/fe- 


\q k ~ 

-5 


11 + 

-5/fc-6 ^ 


: (0 1) 


fl 2) (2 5) (5 6) 


1 r fe i 

+ K}fc- 


r k- 

-1 {9n- 


-i/fe- 


r fc- 

-3 K- 


-3 1 
-3/fc- 


r k- 

-4 t5n- 


-4\ 

-4/ fc _ 6 + 


: (0, 1) 


(1, 3) (3, 4) (4, 6) 


+ {9n} k - 


r fe- 
-1 \9n- 




r fc- 
-3 \Sn- 


-3/fc- 


r fe- 
-5 \9n- 


-5 } fc-6 + 


: (0- 1) 


(1,3) (3, 5) (5, 6) 


+ {9n} k - 


1 k- 
-1 \9n- 


-i}fe- 


r k- 

-4 l-9n- 




1 k- 
-5 \9n- 


-5}fc- 6 + 


i (0, 1) 


(1,4) (4, 5) (5, 6) 


+ {9 h n) k - 


r fe- 
-2 \9n- 


: 2 2 } fe - 


r fc- 

-3 \9n- 


-3/fe- 


1 k- 
-4 \9n- 


-i)k-6 + 


; (0, 2) 


(2, 3) (3, 4) (4, 6) 


+ {9 k n) k - 


r fc- 
-2 \9n- 


-2}*;- 


1 k- 

-3 \9n- 


-3/fc- 


1 k- 
.5 \9n- 


-5}fe- 6 + 


: (0, 2) 


(2, 3) (3, 5) (5, 6) 


+ {s»}fe- 


r fe- 

-2 (Sn- 




-4 i^n- 




r fe- 

-5 \9n- 


-5}fc- 6 + 


(0, 2) 


(2, 4) (4, 5) (5, 6) 


+ {**}*- 


-3 {9n- 


-3JV 


-4 R- 


-4JV 


- 5 


-5/ fc-6 


: (0, 3) 


(3, 4) (4, 5) (5, 6) 



One can construct the sequences of indices (//, v) by considering a line with £ + 1 
points labeled in ascending order from to i. By further considering p unidirectional 
non-stationary jumps from the point labeled to the point labeled £, one gets a 
sequence of initial and final points for the jumps e.g. {(0, 1) (1, 2) (2, 4) (4, 6)} for 
I = 6 and p = 4. We define the ensemble S^ +1 of all such possible jumps, and by 
{Sp +1 }. its i:th element. We can hence write the sum of all terms with p factors 



m 



as 



(3.11) 

whereupon we find 
(3-12) {l k Y k - 



£ n teu 

i (M,-)e{5^+ 1 } i 



-m — 1 



P =i 



With this, we can simply write the now first order recursions eqs. p.lOp as 

(3.i3) mt) - wjh+ew^-^Jh- 



cr=0 
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or using the properties of {7^} 



k—m—X 5 



(3-14) = {fl*}^ + E 



0=1 



! a k -' T - 1 \ 

\i)n-<T-XSk—j ' 



The elimination process can be iterated until only three terms remain in 

n-1 

. -, n— 4 



(3.i5) Uk = £ { 7 ^;::;^ + {r^}"- 4 Mfc 



J=n— 3 

, n— 4 



r t-in— 4 r 1,-, n— 3 r^l-l 71- 4 



-3 

Here {7^}^_^_ 3 , {7n}fe„, 3 l _2 anc ^ {^n}" 4 can ^ e calculated explicitly with the 
help of I gn-u \ ■ Since I 9 n Z^ \ are also known explicitly, eq. (|3.15p is an 

I J k—v I J fc— v 

explicit solution for u k knowing Ufc_„_3, Uk~ n ~2 and for some relevant n. 

The only assumption we have made here is that bi ^ 0, Vi. Unfortunately eq. 
(|3.1o[) is quite cumbersome. It is hence difficult to make any general statement on 
the convergence of eq. (|3.15p . By the simplicity of this scheme, we see no apparent 
reason why it should not be possible to generalize to any finite order recursion 

k 

(3.16) "< »' i 



.jUi-j. 



3.3. An explicit form for A(0). Using the results of section [321 e Q- <|3 -3|> can 
be written 

(3.17) A fe = { 7 *} 



' 1 " ' A fc _„_ 3 + {lnTk-l-2 A fe-»-2 + { r n}" ^ A k-n 



k — n—3 



with 



<T=0 

r<r+i 



(T+l 

p=i 



r fc— a— \\ 
\fln-o— l/fe-j ' 



Using the properties of the coefficients of eq. I|3.3p . 



(3.19){^ +1 },_. = 



-3 
\ n— 5 



n 



0k_(„-3)-0k_l 'llj'=3 



, 1 < J < n — 3 
, j = n-2 

,3 = n-l 



n-3 
—3 ^fe— i 



(3.20) = (-1) 



t»-5 "fc • • a fc-2 



n-3 „2 



fc-1 



n 

J=3 



or 



fc-j-i 



We have thus found an explicit expression for A (0) . Note that the denominators 
diverge for those values of a for which /3j = 1 — a, = 0, Vi € Z. 

n 



Unfortunately the form of eq. Q3.17P does not, a priori, shed light on specific 
properties of the solutions of the Whittaker-Hill formula eq. (|2.13[1 . 

3.4. Numerical considerations. Considering the minors of determinants equiv- 
alent to A„, R. Sips was able to produce a method for calculating the characteristic 
exponent |13| . His method effectively produces first order linear recurrences for n 
minors Di with coefficients an, whereupon the determinant A„ can be written 

A„ = Dl-a\Dl 

For small values of q, his method converges rapidly since on = decreases 
rapidly with growing index i. Having calculated a A„ one has, however, to calculate 
A n / of higher order, to repeat the whole procedure. 

We see that the convergence of our recurrence eq. I|3.3p 

A 4 = ftAj-i - a^A^ + aiO-^Ais, 

is also directly linked to the behavior of oti (ft = 1 — ctfj), The cost of further 
calculating determinants of higher orders will, however, be considerably less than 
for the method of R. Sips, since one only has to successively iterate eq. (|3.3[) - 

4. Conclusion 

With the Whittaker-Hill formula we have derived a third order determinantal 
recursion. By the decomposition of the third order recursion to first order recursions 
we can write any relevant term of the recursion explicitly. We believe that this 
scheme can be generalized to any linear recursion of order greater than one. The 
explicit formula can be used to write the sought after determinant explicitly. We 
expect applications to the numerical treatment of the relation between eigenvalues 
and characteristic exponents. 
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